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Abstract. In this paper we study homoclinic tangles formed by transversal in- 
tersections of the stable and the unstable manifold of a non-resonant, dissipative 
homoclinic saddle point in periodically perturbed second order equations. We prove 
that the dynamics of these homoclinic tangles are that of infinitely wrapped horse- 
shoe maps (See Section [TJ^) . Using as a parameter representing the magnitude 
of the perturbations, we prove that (a) there exist infinitely many disjoint open 
intervals of /x, accumulating at fi — 0, such that the entire homoclinic tangle of the 
perturbed equation consists of one single horseshoe of infinitely many symbols, (b) 
there are parameters in between each of these parameter intervals, such that the 
homoclinic tangle contains attracting periodic solutions, and (c) there are also pa- 
rameters in between where the homoclinic tangles admit non-degenerate transversal 
homoclinic tangency of certain dissipative hyperbolic periodic solutions. In particu- 
lar, (c) implies the existence of strange attractors with SRB measures for a positive 
measure set of parameters. 



1. Introduction 

We start with an autonomous second-order ordinary differential equation tliat con- 
tains a non-resonant, dissipative saddle fixed point with a homoclinic solution. This 
autonomous equation is then subjected to time periodic perturbations. In this paper 
we study the dynamics of the homoclinic tangles formed by transversal intersections 
of the stable and the unstable manifold of the perturbed saddle point. These ho- 
moclinic tangles have been one of the major inspirations for the dynamical systems 
theory and a long standing puzzle in the studies of ordinary differential equations in 
modern times. 

A. Description of results. Instead of focusing on the picture of the globally 

induced time-T maps, by which H. Poincare observed an exceedingly complicated 

mess [P] and S. Smale constructed an embedded horseshoe map [S], we compute 

the return maps induced by periodically perturbed equations around the homoclinic 

solution in the extended phase space. It has turned out that the return maps for these 

homoclinic tangles are infinitely wrapped horseshoe maps, the geometric structure of 

which is as follows. Take an annulus A = x I. We represent points in and 

/ by using variables 6 and z respectively. We call the direction of 6 the horizontal 

direction and the direction of z the vertical direction. To form an infinitely wrapped 

horseshoe map, which we denote as JF, we first divide A into two vertical strips, which 

we denote as V and U. : V ^ A is defined on V but not on U. We compress V 

in the vertical direction and stretch it in the horizontal direction, making the image 
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infinitely long towards both ends. Then we fold it and wrap it around the annulus A 
infinitely many times. See Fig. 1. 

We use a small parameter /i to represent the magnitude of the time periodic pertur- 
bations. Denote the return maps obtained from the periodically perturbed equations 
as J^fj, and letQ 

n^ = {{e,z)eV: T^ie, z)^v Vn > o}, = n„>o 

Then represents all solutions that stay close to the unperturbed homoclinic loop 
in forward times; is the set is attracted to, representing all solutions that stay 
close to the unperturbed homoclinic loop in both the forward and the backward times. 
The geometrical and dynamical structures of the homoclinic tangle represented by JF^ 
are manifested in those of fi^ and A^. A^ obviously contain a horseshoe of infinitely 
many symbols for all /i. This horseshoe covers Smale's horseshoe and all its variations. 
It is the one that resides inside all homoclinic tangles. 




Fig. 1 Infinitely wrapped horseshoe maps. 

The structure of and A^ depend sensitively on the location of the folded part of 
Tp(y\ If this part is deep inside of ?7, then the entire homoclinic tangle is reduced 
to one horseshoe of infinitely many symbols. If it is located inside of V ^ then the 
homoclinic tangles are likely to have attracting periodic solutions or sinks and ob- 
servable chaos associated with non-degenerate transversal homoclinic tangency. We 
prove that, as /i 0, the folded part of T^j^iy) moves horizontally towards Q = +oo 
with a roughly constant speed with respect to p = In fi, crossing V and U infinitely 
many times along the way. It then follows, under mild assumptions, that (a) there 
are infinitely many disjoint open intervals of n, accumulating at = 0, such that the 
entire homoclinic tangle consists of one single horseshoe of infinitely many symbols; 
(b) there are other parameters in between these intervals, such that the homoclinic 
tangle contains attracting periodic solutions; and (c) there are also parameters in 
between where the homoclinic tangle admits non-degenerate transversal homoclinic 

""^We caution that V and U depend also on /i. So to be completely rigorous we ought to write 
and C/p instead of V and U. However, for JF^ derived from the periodically perturbed equations, V 
and U vary only slightly as fi varies, and we could practically think them as being independent of fi 
at this stage. 
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tangency. See Theorems [HIS] in Sect. 14.21 for more details. In particular, (c) implies 
the existence of strange attractors with SRB measures for a positive measure set of 
parameters. 

B. Method of study. We use variables {x,y) to represent the phase space of the 
unperturbed equation and let {x,y) = (0,0) be the saddle fixed point. Denote the 
homoclinic solution for (x, y) = (0, 0) as i. We construct a small neighborhood of i 
by taking the union of a small neighborhood of (0, 0) and a small neighborhood D 
around i out of Ui^. See Fig. 2. Let cr^ e t4 fl D be the two line segments depicted 
in Fig. 2, both perpendicular to the homoclinic solution. We use an angular variable 
^ G S*^ to represent the time. 




D 



Fig. 2 Ue, D and a^. 
In the extended phase space (x, y, 6) we denote 

Ue = UeXS\ D = D X ^1 

and let 

S± = a± X S\ 




Fig. 3 ^^ andM. 

Let Af : S"*" — ^ E~ be the maps induced by the solutions on Ue and : S~ — > S"*" 
be the maps induced by the solutions on D. See Fig. 3. We first compute AA and Af 
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separately, then compose M and M. to obtain an explicit formula for the return map 

We follow the steps of |WOj in deriving the return maps. There are, however, 
two main differences between the classical scenario of homoclinic tangles we now 
consider and the ones studied in |W0] . First, the return maps of this paper are only 
partially defined on E~. After following the entire length of the homoclinic loop of 
the unperturbed equation, part of S~ (represented by V in A) would hit S"*" on one 
side of the local stable manifold of the perturbed saddle where they return to S~; 
and the rest (represented by U in A) would hit on the other side where they sneak 
out. See Fig. 4. Second, analytic controls represented by the estimates in [WO] 
would deteriorate as we approach to the transversal intersections of the stable and the 
unstable manifold of the perturbed saddle, potentially devastating the usefulness of 
the formulas obtained for the return maps. Between the two, the second is essentially 
a technical issue we need to overcome. The first is an intrinsic character of these 
homoclinic tangles. 




Fig. 4 Partial returns to S . 

This paper is organized as follows. In Section [2] we introduce the equations of study 
and a set of changes of variables to transform the equations into certain canonical 
forms. In Section [3] we compute the return maps. Proposition 13.31 in Section [3] is the 
main result of this paper. Theorems about homoclinic tangles are then formulated and 
proved in Section H] assuming the forcing function is in the form of smut. Homoclinic 
tangles associated with general forcing functions are studied in Section [5l 

C. Remarks on history. Homoclinic tangles formed by transversal intersections 
of the stable and the unstable manifold of a periodically perturbed homoclinic saddle 
in systems of ordinary differential equations were first observed by H. Poincare [P] 
more than one hundred years ago. His observation is regarded in general as the event 
that gave birth to the modern theory of chaos and dynamical systems. 

There exists a vast literature on periodically perturbed differential equations (see 
for instance, the reference hst of |GHj ). We could put the related studies roughly 
into two categories. The first category contains the ones that attempted to under- 
stand the dynamics of the associated homoclinic tangles and the second contains the 
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ones that attempted to verify the existence of these homochnic tangles in concrete 
systems of differential equations. Among the most influential in the first category 
are the studies of Cartwright and Littlewood [CL] and Levinson [L] on van der Pol's 
equation and the studies of Sitnikov [SitJ and Alekeseev [A] on Sitnikov's motions. 
These studies led eventually to Smale's construction of his horseshoe map [Sj. To the 
best of our knowledge, Smale's horseshoe is essentially all that has been rigorously 
claimed for these homoclinic tangles. In the second category, the most influential is 
the development of Melnikov's method [M], purposed on verifying the existence of 
homoclinic tangles in concrete systems of differential equations. 

Our idea of constructing return maps such as those derived in |WUj and in this 
paper is motivated by a work of Afraimovich and Shilnikov, see also |AHj . With the 
return maps obtained in [WOj and in this paper, we are able to apply many existing 
theories on maps to the studies of homoclinic tangles of differential equations. Among 
the theories directly applied are the Newhouse theory [N], |PT] on homoclinic tan- 
gency; the theory of SRB measures [Si], [Rj, |Boj : the theory of Henon-hke attractors 
[BCj . |MV] . |BY] : and the recent theory of rank one chaos |WYl] - |WY3j based on the 
theory of Benedicks and Carleson on strongly dissipative Henon maps [BCj . These 
derived maps have led us to many new results. 

In Sect. 14.41 we present an overview on various dynamics scenarios newly found 
around periodically perturbed homolcinic solutions. We now know that, for the two 
main scenarios for the time-T maps depicted in Fig. 5, the dynamics of the first is 
that of an infinitely wrapped horseshoe maps and the second is that of a rank one 
maps. 



In this section we first introduce the equations. We then introduce coordinate 
changes to transform these equations into canonical forms, which we will use in Section 
[3] to compute the return maps. 

2.1. Equations of study. Let (x, y) E be the phase variables and t be the time. 
We start with an autonomous system 




2. Equations and Canonical Forms 



(2.1) 



dx 
'dt 



ax + f{x,y), 



dy 
dt 



Py + 9{x,y) 
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where f{x,y), g{x,y) are real-analytic functions defined on an open domain V C 
satisfying /(0,0) = g{0,0) = dj{0,0) = 5,/(0,0) = d,g{0,0) = dyg{0,0) = 0. First 
we assume that (0, 0) is a non-resonant, dissipative saddle point. To be more precise 
we assume 

(HI) (i) there exists di, d2 > so that for all n,m E Z+, 

\na — m(3\ > di{\n\ + \m\)~'^'^; 

(a) < /3 < a. 

(Hl)(i) is a Diophantine non-resonance condition on a and (3. (Hl)(ii) claims that 
the saddle point (0, 0) is dissipative. Let us also assume that the positive x-side of the 
local stable manifold and the positive y-side of the local unstable manifold of (0, 0) 
are included as part of a homoclinic solution, which we denote as 

i={i{t) = (a(t),6(t)) G V, t G R}. 

Let Q(t) : M — s> M be a real analytic function satisfying Q{t) = Q(t + 2tt). To the 
right of equation (12.11) we add forcing terms to form a non-autonomous system 

dx 

— = -ax + f{x, y) + nA{x, y){p + Q{ujt)), 

(2.2) f 
dv 

-77 = Py + 9{x,y) + fiB{x,y){p + Q{ujt)) 
dt 

where A{x,y), B{x,y) are also real analytic on V satisfying y4(0,0) = -8(0,0) = 0, 
d,A{0,0) = dyA{0,0) = d,B{0,0) = 9,5(0,0) = 0. We regard a, (3, f{x,y),g{x,y), 
A{x,y), B{x,y) and Q(t) as been fixed. uj,p,fi are forcing parameters. The range P 
for {u, p, p) is as follows: Let i?^ >> 1 be an arbitrarily picked real number sufficiently 
large, » Rp » R^. Define 

(2.3) P = {(cu, p, /i) G M^ <uj < R^, Rp^ < p< Rp, 0< p< i?^^}. 

This study is exclusively on equation (12. 2p with parameters inside of P. From Rp » 
Rp» Rui we have, for all (a;, p, p) G P, 

pp, puj « 1. 

We can make the magnitude of \pp\ and \puj\ as small as we desire by adjusting Rp. 

We also fix the values of uj and p, leaving p then as the only parameter to vary. 
We study the solutions of equation (12. 2p in the surroundings of the homoclinic loop 
^ in the original phase space (x, y), which we divide into a small neighborhood Ui. of 
(0,0) and a small neighborhood D around £ out of f/i^. See Fig. 2 in Section [TJ In 
Sect. I2.2l we introduce a change of variables to linearize equation (12. 2p on t/^. In Sect. 
12.31 we introduce coordinate changes to transform (12. 2p on D into a canonical form. 
In the rest of this paper r > is reserved for an integer arbitrarily fixed, and we will 
control the C'-norm of the derived return maps in phase variables and parameters. 

Two small scales: p << e << 1 represent two small scales of different magnitude. 
e represents the size of a small neighborhood of (x, y) = (0, 0), where the linearizations 



of Sect. 12.21 is valid. Define 

and let L+, —L^ be the respective times at which the homoclinic solution l{t) enters 
f/ig in the positive and the negative directions. and are related, both deter- 
mined completely by e and £{t). The parameter n << e controls the magnitude of 
the time-periodic perturbation. We make 

Rf, >> » Rp » R^. 

Notation: Quantities that are independent of the phase variables, time and /i are 
regarded as constants and K is used to denote a generic constant, the precise value 
of which is allowed to change from line to line. On occasion, a specific constant is 
used in different places. We use subscripts to denote such constants as Kq, Ki, - ■ ■ . 
We will also make distinctions between constants depend on e and those do not by 
making such dependencies explicit. A constant that depends on e is written as K{e). 
A constant written as K is independent of e. 

2.2. Linearization on Ue- From this point on all functions are regarded as functions 
in phase variables, time t and the parameter fi. Let X, Y be such that 

x = X + P{X, Y) + fiP{X, Y, e- fi) 

y = Y + QiX,Y)+fxQiX,Y,e;fx) 

where P,Q,P,Q as functions of X and Y are real-analytic on |(X, y)| < 26, and 
the values of these functions and their first derivatives with respect to X and Y 
at {X,Y) = (0,0) are all zero. As is explicitly indicated in (12.41) . P and Q are 
independent of 9 and /z. We also assume that 

P{X, Y,9 + 2tx- ii) = P(X, Y, 9- /x), Q(X, Y,9 + 27r; ^i) = Q(X, Y, 9; /i) 

are periodic of period 27r in 9 and they are also real-analytic with respect to 9 and /i 
for all 6* G M and < R^^- Substituting 9 by ut in (12. 4p defines a non- autonomous, 
near identity coordinate transformation from x, y to (X, Y) , which we write explicitly 
as 

x = X + P(X, Y) + /iP(X, r, ut- /i) 
^^■^^ y = Y + Q{X,Y)+ iJiQ{X,Y,ujt-ii). 

We have 

Proposition 2.1. Assume that a and (3 satisfy the Diophantine non-resonance condi- 
tion (Hl)(i). Then there exists a small neighborhood of {0, 0), the size of which are 
completely determined by equation Ii2.1\) and di,d2 in (Hl)(i), such that there exists 
an analytic coordinate transformation in the form of l{2.5\) that transforms equation 
/ fO)) into 

dX dY 



Moreover, the -norms of P,Q,P,Q as functions of X,Y,6, fi are all uniformly 
bounded from above by a constant K that is independent of both e and fi on {X, Y) G 
Ue, 9eR and fie {-R-\R-^). 

Proof: This is a standard linearization result. See for instance |CLSj for a proof. □ 

2.3. A canonical form around homoclinic loop. In this subsection we derive a 
standard form for equation (12. 2p around the homoclinic loop of equation (12.1 p outside 
oilAi^. Let 

a = {i{t) = (a(t), b{t)) G R^ te M} 
be the homoclinic solution of the unperturbed equation (12.11) . and 



{u{t),vit)) 



be the unit tangent vector of I at ^{t). Let us regard t in [{t) = {a{t),b{t)) not as 
time, but as a parameter that parameterize the curve i in (x, ?/)-space. We replace t 
by s to write this homoclinic loop as i{s) = (a(s), b{s)). We have 

(2.6) ^ = -aais) + f{a{s), bis)), ^ = pb{s) + <?(«(.), bis)). 



By definition, 

u{s) 

(2.7) 

v{s) 



-aa{s) + f{a(s), b{s)) 



^(^-aa{s) + f{a{s), b{s)))^ + {f3b{s) + g{a{s), b{s)))^ ' 
(3b{s)+g{a{s),b{s)) 



^(^-aa{s) + f{a{s), 6(s)))2 + {f3b{s) + g{a{s), bis)))' ' 
Let 

e(s) = {v{s),-u{s)). 
We now introduce new variables {s, z) such that 

{x,y) = i{s) + ze{s). 

This is to say that 

(2.8) X = x{s, z) := a(s) + v{s)z, y = y{s, z) := b{s) — u{s)z. 

We derive the equations for (12.21) in new variables {s,z) defined through (12. 8p . 
Differentiating (12.81) we obtain 

^ = i-aa{s) + /(a(s), bis)) + v'is)z)^ + v{s)^ 

^ = ms) + 9{a{s), b{s)) - u{s)z)^^ - u{s)^^ 
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where u'{s) = ^^^j^, v'{s) = Let us denote 

F{s, z) = — a(a(s) + zv{s)) + f{a{s) + zv{s), b{s) — zu{s)), 
G{s, z) = P{b{s) — zu{s)) + g{a{s) + zv{s), b{s) — zu{s)), 
A{s,z) = A{x{s,z),y{s,z)), 
M{s,z) = B{x{s,z),y{s,z)). 

By using equation fl2.2p . we obtain from equation (12. 9p the new equations for s, z as 
dz 

— = v{s)F{s, z) - m(s)G(s, z) + ^(t;(s)A(s, z) - m(s)B(s, z))(p + Q{yjt)) 

ds _ v{s)G{s, z) + u(s)F(s, z) + /i(t;(s)B(s, z) + m(s)A(s, z)){p + Q(a;t)) 
c^t ~ ^F{s, 0)2 + ^(s, 0)2 + z{u{s)v'{s) - v{s)u'{s)) 

We re-write these equations as 



(2.10) 

where 
(2.11) 



dz 

— = E{s)z + z^W2{s, z) + fi{v{s)A{s, z) - u{s)M{s, z)){p + Q{ujt)) 

ds fi{v{s)M{s,z)+u{s)A{s,z)){p+Q{u;t)) 
1 + zwi[s, z, out; jj.) + 



dt ' ' ' ' y/Fis,Oy + G{s,0) 
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Eis) = v'is)i-a + dJiais), b{s))) + u\s){P + dygiais), 6(s))) 
- u{s)v{s){dyf{a{s), b{s)) + d^g{a{s), b{s))). 



and Wi{s, z,9 + 2n; jj) = Wi{s, z, 9; fj,) is periodic in 9 of period 27r. Also in the rest of 
this section we let Ki{e) be a given constant independent of /i, and regard equation 
(I2.10p as been defined on 

{s G [-2L-, 2L+], 1^1 < K^{e)fi; /i G (0, R-^')}. 

The C""-norms of Wi{s, z, 9; jj) and ^2(5, z) are bounded above by a constant K{e). 
Finally we re-scale the variable z by letting 

(2.12) Z = ii-^z. 

We arrive at the following equations 
dZ 

— = E{s)Z + fiW2{s, Z, ut; fi) + H(s)(p + Q{ujt)) 

(2.13) ^1 

— = 1 + /iWi(s, Z,Ljt;fi) 
dt 

where 

(2.14) H(s) = v{s)A{a{s), b{s)) - u{s)B{a{s), 6(s)); 

and (s, Z] /i) are defined on 

© = {(s,Z;/i): sG [-2L-,2L+], \Z\ < K,{e), e {0, R~')}. 
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Note that here we assume that R^^ is sufficiently small so that 

/i<< min (F(s,Of + G(s,OY). 

se[-2L- ,2L+] 

Again, the norms of Wi, W2 with respect to s, Z; ^ and t are uniformly bounded 
by a constant K{e) for (s, Z; /i) G D and t G M. Equation fl2.13p is the one we need. 

2.4. Technical estimates. Estimates presented in this subsection are directly taken 
from [WO], which we include for completeness. 

Notation: We are going to adopt the following convention in comparing the mag- 
nitude of two functions f{t) and g{t). We denote f{t) -< g{t) if there exists > 
independent of t so that \ f{t)\ < K\g(t)\ as t ^ oo (or — oo). We denote f{t) ~ g{t) 
if in addition we have > K^^\g{t)\. We also denote f{t) ^ g{t) if 

m . 

git) 

as t ^ oo (or — oo). 

Recall that l{t) = {a{t),b{t)) is the homoclinic solution for the hyperbolic fixed 
point (0, 0) of equation (12. ip . {u{t),v{t)) is the unit tangent vector of i at £{t) defined 
through fl^ . 

Lemma 2.1. As t —* +oo, 

a{t) ~ e""*, b{t) ^ e-2°*, u{t) ^ -1, v{t) ^ e""*; 
a{-t) -< e~^^\ b{-t) ~ e-"', u{-t) -< e"^*, v{~t) ^ 1. 

Proof: We are simply re-stating the fact that i{t) (0, 0) with an exponential rate 
—a in the positive time direction along the x-axis, and an exponential rate f3 in the 
negative time direction along the y-axis. □ 
Let E{s) be as in fl2TTD . 

Lemma 2.2. As +oo, 

(i) f_^_{E{s) + a)ds -< 1, jl;^ {E{s) - (3)ds -< 1. 

(ii) /°^_ E{s)ds ^ ~aL-, J^'' E{s)ds ^ (3L+. 

Proof: (i) claims that the integrals are convergent as oo. For the ffist integral, 

we observe that by adding a to E{t), we obtain E{t) + a as a collection of terms, 
each of which decays exponentially as t ^ — oo according to Lemma [2.1[ Similarly, 
taking (3 away from E{t), we obtain E{t) — /3 as a collection of terms, each of which 
decays exponentially as t — >• oo. 
For (ii) we write 

E{s)ds = -aL~ + / {E{s) + a)ds 

L- J-L- 

E{s)ds = (3L++ / {E{s)-f3)ds. 
'0 Jo 
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(ii) now follows from (i). □ 
We also have 

Lemma 2.3. As e ^0, e e~°'^^ ~ e"^^". 

Proof: This follows directly from the definition of and Lemma 12. 1[ □ 

3. Derivation of Return Maps 
Let 6 E he an angular variable for time. We re-write equation (12.21) as 
dx 

— = -ax + f{x, y) + y){p + 2(6*)), 
at 

(3.1) § = + + y^^p + 

de 

— = OJ. 

dt 

S^*" are formally defined in Sect. 13.11 (See Section [TB). In Sect. 13.21 we study coordi- 
nate conversions between (X, F, 9) and (s, Z, 9) on S^. M : S"*" — > S~, : E~ — »• E"*" 
and the return map J-' = N o M. : are computed in Sect. 13. 3[ 

3.1. Poincare sections S^. We start with equation 03.11) on U^. We have obtained 
in Sect. 12.21 a change of variables on in the form of (12.41) . that is, 

x = X + P(X, F) + /xP(X, F, ^; /i) 
l/ = F + g(X,F)+/iQ(X,y,^;/x) 

that transforms equation (13.11) to the linear equation 
/ X dX dY d9 

on We define S^*" inside of fl D by letting 

s- = {(X, r, 0) : r = £, |x| <ij, 9e s^}, 

and 

S+ = {(X, F, ^) : X = e, |F| < iri(£)/i, 9 G S^}. 

-ft'i(£) will be precisely defined in Sect. 13.31 Observe that in [WO], S^*" are defined 
in slightly different terms. The current definition is designed to avoid the long and 
deteriorating derivative estimates of [WO]. 

We turn to the canonical form for equation (13.11) on D. Let 

(3.4) X = a{s) + fiv{s)Z, y = b{s) — fiu{s)Z. 
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Then according to Sect. I2.3[ equation (13 .ip on D is written in {s, Z,9) as 

— = E{s)Z + fiW2{s, Z, 6- /i) + H(s)(p + Q{9)) 
at 

(3.5) ^ = l + /xu;i(s,Z,^;/i) 

dt 

dO 

— = ijj 
dt 



where 
(3.6) 



E{s) = v\s){-a + dJiais), b{s))) + u^{s){P + dyg{a{s), bis))) 
- u{s)v{s){dyf{a{s), b{s)) + d^g{a{s), b{s)))] 



(3.7) H(s) = v{s)A{a{s), b{s)) - u{s)B{a{s), b{s))] 

and the C^-norms of Wi, W2 are bounded from above by a K{e) on D x (0, R^^) where 
D = {{s,Z,9) ■ s e [-2L-,2L+l \Z\ < Ki{e), 6 e S^}. 

Let q G S+ or S^. We represent q by using the (X, F, 6*) -coordinates, for which 
we have X = e on S"*" and Y = e on S~. We can also use (s, Z, 0)-coordinate to 
represent the same g, for which the defining equations for Yj^ are not as direct. To 
compute the return maps, we need to first attend two issues that are technical in 
nature. First, we need to derive the defining equations for Yj^ for {s,Z,9). Second, 
we need to be able to change coordinates from (X, Y, 9) to (s, Z, 9) and vice versa on 
S^. We start with some preparations in notation. 

Notation: The intended formula for the return maps would inevitably contain 
terms that are explicit and terms that are implicit. Implicit terms are usually "error" 
terms, and the usefulness of a derived formula would depend completely on how 
well the error terms are controlled. In this paper we aim on C""-control on all error 
terms. The derivations of the return maps involve a composition of maps and multiple 
coordinate changes. To facilitate our presentation, from this point on we adopt specific 
conventions for indicating controls on magnitude. For a given constant, we write 0{1), 
0{e) or 0(yu) to indicate that the magnitude of the constant is bounded by Ke 
or K{e)fi, respectively. For a function of a set V of variables on a specific domain, 
we write Ov{l), Ov{£) or Ovifi') to indicate that the C^-norm of the function on the 
specified domain is bounded by K,Ke or K{e)fi, respectively. We chose to specify 
the domain in the surrounding text rather than explicitly involving it in the notation. 
For example, Oz^ifi) represents a function of Z, 9, the C"'-norm of which is bounded 
above by K{e)fi. 

The new parameter p: For the formulas obtained to be most useful, it is also 
desirable that we have control on the derivatives with respect to the forcing parameter 
IJ,. Taking derivative with respect /i, however, are problematic because such action 
takes fi out of the needed places. To resolve this potentially damaging problem we 
introduce a new parameter p = In /i and regard p, not /x, as our bottom-line parameter. 
In another word, we regard shorthand for e^, and all functions written in fx 
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as functions in p. Observe that /i G (0,yUo) corresponds to j9 G (— cxo, In/io). This 
is a very important conceptual point because by regarding a function -F(yu) of /i as a 
function of p, we have 

^,F{^^)=^^^,F{^^). 

So regarding as a function of p would give back to us that much needed factor 
/i in derivative estimates. 

3.2. Conversion of coordinates on S^. We start with the defining equations for 
S"*" in {s,Z,6). Results for are similar. 

Lemma 3.1. We have for {s, Z, 6) G S"*" 

s = L+ + Cz,e,p(/i). 

Proof: From fl3.2p and fl3.4l) . we have on 

a{s) + v{s)z = e + P{e, Y) + iJ,P{e, Y, 6; /i) 

bis) - u{s)z = Y + Q{e, Y) + fiQ{e, Y, 9; /i). 



(3.8) 
By definition 
(3.9) 
Let 
(3.10) 



a(L+) =e + P{e,0) 
b{L+) = QiE,0). 



Wi = a{s) - a{L^) + v{s)z - ijP{e, 0, 6; jj), 
W2 = b{s) - 6(L+) - u{s)z - fiQ{e, 0, 6; jj). 
We have from (ESD and M . 

VTi = P(£, Y) - P{e, 0) + /i(P(£, r, 0; /i) - P{e, 0, /i)) 
1^2 = 1^ + Q{e, Y) - g(£, 0) + /i((Q(£, y, ^; /i) - Q{e, 0, ^; /i)) 
which we re-write as 

^ ^1 = ((^(^) + /^a,p(l))V' + Oy,e,,m' 

W2 = {1 + 0{e) + fiOeM)Y + OY,eA'^)Y\ 

We first obtain 

(3.12) F = (1 + 0(e) + ^lOe,p{l))W2 + 0^,,,,p(l)l^| 

by inverting the second line in (13.111) . We then substitute into the first line in (13.110 
to obtain 

Wi = {0{e) + /ia,p(l))((l + 0{e) + /ia,p(l))W^2 + Ow„e,p{'^)Wi) 
+ Oy,e,p(l)((l + 0{e) + /ia,p(l))1^2 + Ow^„e,p(l)W^|)' 
= {0{e) + ^a,p(l))W^2 + Ow,Av{^)Wl 
Consequently, 

(3.13) F(s, Z, e, /i) := H^i - {0{e) + ^a,p(l))1^2 + = 0, 
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where Wi, W2 as function of s, Z, 9 and /x are defined by fl3.10p . To re- write Wi, W2 
we let 

(3.14) ^ = s-L-^ 
and expand a(s) in ^ as 

00 

a{s) = a(L+) + a'{L+)^ + ^ a,{L+)C. 

1=2 

Expansions for b{s),u{s) and v{s) are similar. We have 

00 00 

i=2 i=2 

— uP(e, 0, 9; 11) 

(3.15) P V > > 

W2 = h'{L+)i + HL^)C - u{L^)z - {u\L^)i + 5^ «.(L+)e)z 

j=2 i=2 

- fiQ{e,0,9;fi). 

We now put f lS.lSp for 14^1, W2 back into equation f l3.13p and replace z by /zZ. We 
obtain 

(a'(L+) - 0(£)6'(L+) + /i(^^,p,00e = Oz,eAf^) 



where the C" norm of h{9,p,^) is bounded from above by K(e). From Lemma \2.1\ 
a'{L+) ^ -ae, b'{L+) = ole^). We finally obtain 

by solving C,- This proves Lemma [3. II □ 
Lemma [3. II is not precise enough. We need the following refinement. 

Lemma 3.2. We have on 

Proof: It suffices for us to drop all terms that is Oz,e,p{fJ''^) in equation (13.131) to 
solve for ^. From Lemma [3.11 we conclude that all terms in ^, 2; of degree higher than 
one are Oz,e,p{fi'^)- With these terms all dropped, fl3.13p becomes 

(3.16) (a'(L+) - 0{e)b\L+m + (^(^'") + 0{e)u{L+))z = /iO,,p(l), 
from which the estimates of Lemma 13.21 on S"*" follows. □ 
From this point on we let 

X = ij,-^X, Y = ii~^Y. 

Lemma 3.3. On S"*" we have 

Y = (1 + 0{e))Z + OeM + Oz,eAf^). 
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(3.17) 



Proof: We have 

F = (1 + 0{e)){h'{L+)i - u{L^)z - ^iQ{e, 0, 9- /i)) + OzfiA^^^) 

= (1 + 0{e))z + + 

where the first equahty follows from using (]3.12p , ( 13.150 and Lemma I3.lt the second 
equality from using Lemma [3.2[ To obtain the third equality we use u{L~^) = — 1 + 
0{e), a'{L+) ^ -ae, b'{L+) = 0{e^). □ 

Along similar lines we can also prove 

Lemma 3.4. On S^, we have 
(i) s = -L~ + Oz,e,p{f^); and 
(tt) Z = {1 + 0{e))X + Oe,p{l) + O^.eM- 

Proof: Left to the reader as an exercise. □ 

3.3. The return map JF = MoM.. First we compute TV : S+ and : 

S"*" separately. We then compose M and M. by using Lemmas 13.31 and 13.41 

A. The induced map AT : S+ ^ S". For (X, Y, ^) G S+ we have X = e^-^ by 
definition. Similarly, for (X, Y, G S~ we have Y = e/x"^. Denote a point on S+ by 
using (Y, 9) and a point on S~ by using (X, 9), and let 

(Xi,^^i) =N{Y,9) 

for (Y,0) G S+. 

Proposition 3.1. We have for {Y,9) G S+, 

Xi = (/i£-i)t-iYt 

Proof: Let T be the time it takes for the solution of (13. 3p from {e, Y, 9) G S+ to get 

to (Xi,£:,6'i) G We have 

Xi = ee""^, e = ye^"^, = ^ + ujT, 
from which (I3.18P follows. □ 

B. The induced map Al : E~ — S"*". Let EI(s) be as in ( 13. 7p . In what follows, we 
write 



(3.19) 



{9) = J m(s)Qi9 + ujs + tuL-)e~^o'^(^^'^^ds 
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We also write 

(3.20) Pl = e^-^- P+ = e-fo^ ^(^^'^^ 

Note that for Pl we integrate from s = —L~ to s = , while for P^ the integration 
starts from s = 0. First we have 



Lemma 3.5. 



Pr. ^ << 1, Pr"*" ~ £ a >> 1. 



Proof: Both estimates follows directly from Lemmas 12.21 and 12. 3[ □ 

For q = {s~,Z,9) G the value of s~ is uniquely determined by that of {Z,9) 
through Lemma I3.4( i). So it is allowed for us to use {Z,9) to represent q. Let 
(s(t), Z(t), 9(t)) be the solution of equation (13.51) initiated at (s~, Z, 9), and i be the 
time {s{i),Z{i),9{i)) hit S+. By definition M{q) = {s{i), Z{i),9{i)). In what follows 
we write 

s+ = s{i), Z = Z{i), 9 = 9{t). 
Proposition 3.2. Denote {Z,9) =M{Z,9). We have 

Z = PUpAl + <I>l{9)) + PlZ + Oz,eM 

9 = 9 + u;{L-^ + L-) + Oz,eAp). 
Proof: Let us re-write equation (13. 5p as 

-J- = E{s)Z + M{s){p + S(»)) + 0,,z.eM 

(3,22) 

— = u; + C's,z,e,p(yu) 
as 

on D X (0, R~^) where 

D = {(s, Z, 0) : s G [-2L-, 2L+], \Z\ < Ki{e), 9 G S^}. 
Dropping all error terms in (I3.22p we have 

'^^ = E(a)Z + ns)(p+Q(0)) 
(3.23) ^^^^ 

ds 

We estimate the solution of equation (13.22^ initiated at (Z, 9) from s = to s = 
by the solution of equation (I3.23P initiated at the same (Z, 9) from s = —L~ to 
s = L+. By the smooth dependencies of solutions with respect to equations and 
initial conditions, the error of such estimates, according to Lemma 13.11 and Lemma 
Eaii), is 

provided that both solutions stay inside of D. By solving (13.231) . we obtain 
Z = Pl{Z + <^Lm + Oz,eM + O^^eM 

(3.24) 

9 = 9 + uj{L+ + L-) + OzfiM + Oz,eM 
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where Pl is as in (13.201) and 



(3.25) <^^{e) = j H(s)(p+ Q(0 + cuL- +cjr)) ■e-^-V^(^)'i^c/r. 

From (I3.24P we have 

e = e + uj{L+ + L~) + Oz,eM. 



(3.26) 



Let 

(3.27) K,{e)= max P,(2 + |$,(^)|) 

6651, se[-2L-,2L+],/^e(-ij;:\ij;:^) 

where Pg and $<( are obtained by replacing with s in Pl and $l. Ki{e) is the one 
we use for D and S"*". Solutions of (I3.22p initiated on E~ will stay inside of D before 
hitting S^. To finish, we observe that 

= P+(pA^ + 0z.W). 

This finishes the proof of Proposition 13.21 □ 

C. The return map = M o M. We are now ready to compute the return map 
= M o M. : — > S^. We use (X, 6') to represent a point on S~ and denote 
(X,^) = J^(X,^). 

Proposition 3.3. The map = M o M. : — > is g'zven 6y 

X = (^£-i)f-i[(l + 0(£))P+F(X,^)]f 
^^•^^^ ^ = + cu(L+ + L-) + I In /x-^£(l + 0(£))P+ - ^ lnF(X, 6) + Ox,e,p(/u) 

where 

F(X, e) = {pAl + + + 0{e))^ 

+ {Pt)-\l + PL)Oe,,{l) + 0^,eM, 



and Pl, P^ and 0l(6') are as in Ii3.iy\) and ^ 



Proof: By using Proposition 13.21 and Lemma [3.41 we have 

Z = Pl{1 + 0{6))X + P+{pAl + MG)) + PLOeM + O-^^eM 
e = e + iu{L+ + L-) + Ox,eAf^). 
Let Y be the Y-coordinate for {Z, 9), we have from Lemma [3.31 

Y = (1 + C(£))P+F(X,^) 
where F(X, 6) is as in (13:291) . We then obtain ([3:28|) by using (13181) . □ 
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We remark that = M o M. is only defined on the part of S where 

F(X,0) > 0, 

and the set in E~ defined by F = is on the stable manifold of the saddle (x, y) = 
(0, 0). Proposition 13.31 is the main result of this paper. 

4. Dynamics of Homoclinic Tangles: Q{t) = sint 

In this section we let Q{t) = sint in equation (12. 2p . In Sect. I4.1l we derive the return 
maps. In Sect. 14. 2[ we prove that these return maps are infinitely wrapped horseshoe 
maps (See Section [ijA.)- In particular, we prove that (a) there exist infinitely many 
disjoint open intervals of fi, accumulating at = 0, such that the entire homoclinic 
tangle is one single horseshoe represented by a full shift of infinitely many symbols 
(Theorem!!]); (b) there are parameters in between each of these intervals, such that the 
homoclinic tangle contains attracting periodic solutions (Theorem [2]); and (c) there 
are also parameters in between where the homoclinic tangle admits non-degenerate 
transversal homoclinic tangency (Theorem[3]). The existence of Henon-like attractors, 
following directly from Theorem [3] and |MVj . is stated in Corollary 14.11 In Sect. 
14.31 we study the associated homoclinic tangles by numerically iterating the derived 
return maps. Finally in Sect. 14.41 we summarize various dynamics scenarios in the 
surroundings of periodically perturbed homoclinic solutions newly found through the 
return maps of Proposition 13. 3[ 

4.1. The return maps for homoclinic tangle. Let Q{t) = sint in equation (12. 2p . 
Let F(X, 6) be as in Proposition 13. 3[ The stable and the unstable manifold of (x, y) = 
(0, 0) of equation (12.21) intersect if and only if there exists 9 such that F(0, 9) = 0. In 
[WO], the authors excluded the possibility of these intersections by restricting to a 
specific range of forcing parameters. We now allow F(0, 9) = 0. 
Let 




Recall that i{s) = (a(s), b{s)), s G M is the homoclinic solution of equation (12.11) and 
(m(s), f (s)) is the unit tangent vector of £{s). Also recall that M{s) is as in (12.141) and 
E{s) is as in (12.111) . Using the conclusions of Sect. 13.21 it is easy to verify that A, C 
and S are all well-defined. In the rest of this section we assume that 

(H2) (i) Ay^O] and (ii) C\uj) + S\u) ^ 0. 

For a given equation (12.11) satisfying (HI), (H2)(i) holds for majority of A[x,y) 
and B{x,y). (H2)(ii) requires that, as a function of s, the Fourier spectrum of the 
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function 

is not identically zero. We know that R{s) decays exponentially as a function of s, 
and it follows that the Fourier transform R{u)) is analytic in a strip contain the real 
tu-axis. Consequently, R^uj) = for at most a discrete set of values of u because R{s) 
is not identically zero. Note that R{u) = C{uj) + iS{uj). 

Specifications on parameters: The parameters uj,p,e,fi are specified as follows. 
First we fix (arbitrarily) an u such that (H2)(ii) holds. Then we fix a value of p such 
thatS 

3 < -\jC\uj) + S^(uj) < 9. 
pA 

Numbers 3 and 9 here have no specific meaning and can be replaced by any other 
two numbers larger than 1. We let e be small enough for a variety of reasons: one is 
to validate the derivations of the previous sections and another is to make 



(4.2) 2<-^^^Cl{u;) + Sl{u;)<10 

where Al,Cl,Sl are obtained by replacing the integral bounds ±oo with iL^*" re- 
spectively in A, C, S. p {<< e) is the only parameter we allow to vary. 

The return maps: In the rest of this section we use z for X, A for S^. So we write 

A = {{e, z): d e R/(27rZ), \z\ < 1}. 

We regard u, p, e as been fixed. Let (6*1, Zi) = J-'{9, z) for (6*, z) E A where T is from 
Proposition 13. 3[ We have 

e^ = e + ai- %hiW{d,z,p) 

(4.3) P 

z^ = h[w{e,z,p)fp 

where 

a = In/i-i + ^(L+ + L-) + ^ ln(£(l + 0{e))PtALp) 

(4.4) P P 

h={ps-^)^-\{l + 0{e))PtALp\'l 

and 

(4.5) F(^,2,/i) = l + csin^ + k;z + E(^,/i) + Oe,,,p(/i), 
in which 

(4.6) 



{A,p)-\Cl + S^ 



L 

\.= {A,p)-'P,{Pi)-\l + 0{e)) 
and 

(4.7) E(e,/.) = {ALp)-\Ptr\l + PL)Oe,,{l). 

^Let us assume ^ > here to maintain a positive range for p. If A < we need to switch p to 
-P- 
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Note that in getting (14.31) we have changed 6 + ujL + cq to 6 where Cq is such that 
tan Co = Cj^^Sl. a, b, c, k and E(^, /i) follows: 

(i) b ^ as /i ^ 0. We can think JF as an unfolding of the ID maps 

f{e) = ^ + a - ^ ln(l + c sin ^ + E(^, 0)). 

(ii) a +00 as — > 0. a is a large number. But since it appears in the angular 
component we can module it by 27r. With u, p and e been fixed, a is essentially 
a;/3~^ ln/i~^. Varying p from a small /xq > to zero is to run a over (ao,+oo) for 
some ao ~ In/Xg ^ 

(iii) By (14.21) . c G [2, 10] is a constant independent of /i. Consequently, there exists 
an interval for 6 so that 1 + c sin 6 < 0, and the stable and the unstable manifold of 
the perturbed saddle of equation (12.21) do intersect. Also observe that from Lemma 
13.51 we have 

When e is sufficiently small, E(^, /i) is a C'-small perturbation to 1 + c sin 6. 

(iv) k is a small number independent of fi. In fact, k ~ e"^ ^ from Lemma [3.51 k 
is, however, much larger than /i and it follows that the first derivative of F(^, z, p) 
with respect to 2; is ~ k. This implies that the unfolding from f{6) in (i) to T is 
non-degenerate in 2;-direction, and is controlled completely by the linear term kz. 

(v) It is important that E(^, /i) is independent of z. Otherwise we would have 
trouble in controlling what happens in z-direction. See (iv) above. 

New notation on parameter: In the rest of this section we put a in the place of p, 
regarding it as the bottom line parameter. Both /i and p are regarded as functions of 
a. Since we have fixed a;, p and e, c and k are fixed constants independent of a. We 
denote the return maps as Ta to emphasize that a is the parameter, b is a function 
of a. Because a and p are hnearly related, we have Oe^piX) = Oe,a{X) (14.71) and 
Oe,zM = Oe,,M in (BSD . 

4.2. Homoclinic tangles as an infinitely wrapped horseshoe map. For q = 

{9, z) G A, let V = (m, v) be a tangent vector of ^ at g and let s(v) = vu~^. s(y) is the 
slope of V. We say that v is horizontal if |s(v)| < and v is vertical if |s(v)| > 100. 
A curve in .4 is a horizontal curve if all its tangent vectors are horizontal and it is a 
vertical curve if all its tangent vectors are vertical. A vertical curve is fully extended 
if it reaches both boundaries of A in ^-direction. A region in A that is bounded by 
two non-intersecting, fully extended vertical curves is a vertical strip. For a given 
vertical strip V , a horizontal strip in is a region bounded by two non-intersecting 
horizontal curves traversing V in 6'-direction. 
Observe that 

(4.8) F(0, 2, /i) = k^ + 1 + c sin ^ + E(0, p) + Oe,zM = 

defines two fully extended vertical curves that divide A into two vertical strips, which 
we denote as V and U. Let F > on \^ and F < on U. J-g is well-defined on V but 
not on U. U is the window through which the solutions of equation (12. 2p sneak out. 
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Let 



(4.9) n^ = {{e,z)ev: j=i{Q, z) G Vn > 0}, A, = ^„>o^:(^^a). 

f2a represents all solutions of equation (12.21) that stay close to I in forward times; Aa 
is the set fia is attracted to, representing all solutions that stay close to £ in both 
the forward and the backward times, fia and Aa together represent the homoclinic 
tangles, the structure of which we now unravel through J-g. 
For a fixed z G [—1, 1], let 

= {0 G (-^TT, ^vr] : {Q.z)^V}. 

is an interval in (— |vr, |7r), which we denote as {Oi{z), Or{z)). Let = {{9, z) : 9 & 
h}- ^aihz) is a ID curve in A parameterized in 6*, which we denote as (^1(6*), 6*1 (6*)). 
By definition 

(4.10) 0i(0) = + ^lnF(^,z,/i). 
We have 

Lemma 4.1. Assume that uj(3~^ > 100. 

(a) \img^g^(^^-)-{9i,zi) = \img_g^(^^)+{9i, zi) = (+00,0). 

(b) For every fixed z G [—1, 1], there exists a unique value of 9, which we denote as 
9c{z), such that 



(c) Let 

(4.11) Vf = U,e[-i,i]mz)eV 



d9 

d9i 



<2}. 



d9 , 

Then Vj is a vertical strip, the horizontal size of which is < lQuj^^(3. 

Proof: Observe, from (14.81) . that 9i G (— |7r,0) where cos6' > 0, and 9r G (tt, |7r) 
where cos 6' < 0. (a) follows directly from the fact that, as ^ 9^,9~, F ^ 0. To 
prove (b) we first observe that, because F ^ as 6^ ^ 6*^, 



11 + csinr| < JsTe^""' << 1 



and it follows that 
Consequently, 

Similarly, we have 



d¥ 



C COS 9r < — 1 . 



V d9i / ^-il-^F^ 

hm — — = hm 1 — 00 p —t^t; = +00. 



d9 0^0- \ ' ¥89 



hm -777= hm { 1 - up -— = -00. 
■ d9 e->e+ \ ¥ d9 ' 



"I " "I 



21 



Therefore there exists at least one 9c{z) satisfying ^ = 0. For the uniqueness we 
observe that 

rF- h« + csine > 



^2 \QQ2 \de ¥■ 



for all 9. Recall that c > 2. 

To prove (c) we observe that the boundary of Vf is defined by 

from which we obtain 

Icos^l < -u-^(3 + Ke'^''~\ 
2 

(c) follows directly from this estimate. □ 
We are now ready to formally state and prove the first of our theorems. 

Theorem 1 (Horseshoe of infinitely many symbols). Let Q{t) = sint and 
assume (HI) and (H2) for equation h2.2) . Let the parameters uj,p,e been specified as 
in Sect. \4.1\ If in addition ujf3~^ > 100, then there exists a sequence of fi, which we 
denote as 

1 » > /if^ > • • ■ > fx^''^ > fx^J^^ > ■■■>0 
such that for all jj, G [iin\ jj^n'], on 

K = {{e,z)eV ■ Ti{e, z) G l^, Vi G Z} 
conjugates to a full shift of countably many symbols. 

Proof: For different values of /i, the corresponding vertical curves in A defined 
by fl4.8p are 0{fi) close. So V and U are almost stationary as a varies from ag to 
+00. On the other hand, it follows from (14. 3 P that, by varying a from ao to +oo, 
we move J-'aiV) horizontally towards 9 = +oo. Denote T = and let Vf be the 
vertical strip defined through (14. lip . The horizontal size of J^iyf) is smaller than 
20/3u;~^ from Lemma 14.11 assuming > 100, which is in turn smaller than the 

horizontal size of U. Therefore J^{Vf) traverses A infinitely many times in horizontal 
direction as we vary a from ao to +oo and there are infinitely many sub-intervals of 
a, such that J^{Vf) C U. For these parameter values J-'iV) fl V consists of countably 
many horizontal strips in V (see Fig. 1 in Section [TJ/\.), to each of which we assign a 
positive integer according naturally to the order in which these strips are stacked in 
the downward ^-direction. 

For q & A, let v be a tangent vector at q. Let Ch{q) be the collection of all v 
satisfying |s(v)| < ji^, and C„(g) be the collection of all v satisfying |s(v)| > 100. To 
prove that A conjugates to a full shift of all positive integers, it suffices to verify that 
we have, assuming J-'iVf) C U, 

(i) DJ^iChiq)) C ChiJ'iq)) on T-^J'iV) H V), and 

(ii) DJ^'\C,{q)) C C^iJ'iq)) on J^{V) H V. 
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To prove (i) we first compute DJ^ by using (14. 3p . Let (6'i, Zi) = J-'{9, z), we have 

_ -1 1 OF 

(4.12) 



86>i 86>i 
ozi ozi 

ae dz 



1 ,,/Q-119F 



, 1 aF 



dz 



CCOS^^ + £^"" a,a(l) + a,.,a(/i) 



where F = ¥{6, z, /i) is as in (14. 5 p and 

09 

d¥ 

— = k + Oe,^,M- 
Let V be such that |s(v)| < we have from fl4.12p 



(4.13) 



|.(D^(v))| 



lldF, 



We have two cases to consider. 

Case 1: ¥ > -\/k- In this case we have 

From {6, z) G T-^{Tiy) n l^) and J^(Vf) C f/, it follows that {Q, z) ^ Vj therefore 

5-1 1 d¥ 







de 





¥de 



> 2. 



These two estimates together implies that the denominator for |s(DjF(v))| in (14.131) 
is > 1, and it follows that |s(DjF(v))| < 

Case 2: ¥ < Vk.. In this case 

|1 + csin6'| < Ke^ + Vk, 

from which we have 

(4.14) |ccos^|>l. 

It then follows that the denominator for |s(DjF(v))| in (14.131) is > which implies 
|s(D^(v))| < jig. This finishes our proof for (i). 
To prove (ii) we let v be such that |s(v)| > 100. From (14.121) . 

1 d¥ 

(4.15) DJ^-^ 



1 



ap-^h¥^''-'-'f 



and we have 



\s{DT-\m 



-«r^bF"^'^-if .-i(v) + (1 - o^r^if 



a/3-ibF"/3-'-if:s-i(v) - 



F dz 
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We again divide into the cases of F > and F < v^. For the case of F > v^, the 
magnitude of the denominator << 1 and that of the numerator is > 1 again because 



1 d¥ 



> 2 



from the assumption that ( 
as 

\s{DJ^-'M)\ 



^ Vf. For the case ofF < Vk, we re- write \s{DJ^~^{v))\ 



The denominator is again << 1 and the dominating term in the numerator is 



> 1. 



The last estimate is from 



d¥ 



de 



|ccos6'| > 1 



again by fl4.14p . This proves (ii). □ 

We refer the reader to Chapter III.l of [Mo] for a detailed discussion on horseshoes 
of infinitely many symbols. 

Remarks: 1. For the parameters of Theorem [1], the entire homoclinic tangle 
consists of a single horseshoe of infinitely many symbols. 

2. Due to the expansions associated with singularities of the logarithmal function 
in fl4.3p . jFjj induces a horseshoe of infinitely many symbols on V\Vf for all < fiQ. 
This horseshoe covers Smale's horseshoe and all its variations. It is the one that 
resides inside all homoclinic tangles. 

3. A is much more comphcated when J-'a.iVf) intersects V. As J-'aiVf) traverses 
V, we encounter complicated dynamical patterns caused by our allowing the images 
of the unstable manifold of the horseshoe in V \ Vf (see remark 2) to come back to 
traverse the stable manifold of the same horseshoe. We will prove, momentarily, that 
there are parameters that admit periodic sinks and there are also others that admit 
non-degenerate transversal homoclinic tangency. We unfortunately do not have a 
bifurcation diagram for JF^,. However, we know from fl4.3p that the same diagram are 
repeated infinitely many times a.s fi —>■ 0. 

4. We caution that, though the horseshoe of Theorem [1] represents all solutions 
of the perturbed equation that stay forever inside of a small neighborhood of the 
homoclinic loop i, solutions sneaked out through U might find a way to come back to 
A, creating more complicated structures. One particular mechanism for such coming 
back is for the unperturbed equation to have two homoclinic solutions. See Fig. 6(a). 
In this case, part of U would come back to A following the other homoclinic loop. On 
the other hand, it is easy to obtain examples for which the solutions sneaked out of 
U would never come back. In this case the entire homoclinic tangle for the perturbed 
equation is in fact reduced to the horseshoe of Theorem [TJ all it takes for this to 
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happen is for us to send the other branch of the local unstable manifold of (0, 0) to a 
sink. See Fig. 6(b). 





(a) (b) 

Fig. 6 (a) U gets back to A, and (b) all points in U approach a sink. 

Our next Theorem is about the existence of periodic sinks. We remark that these 
periodic sinks are not Newhouse sinks associated with homoclinic tangency. 

Theorem 2 (Periodic sinks). Let the assumptions be identical to that of Theorem 

U\ Then there exists an open set of fi inside each of the intervals [fin \ f^n+i], such 
that the corresponding homoclinic tangle admits a periodic sink. 

Proof: Let 9c{z) be as in Lemma [4.1( b). To make the dependency on /i explicit we 

write it as 6c{z,fi). Let be the value of a at ;U = fin^ and [a„] = a„ — a„mo(i(27r). 

Observe that there exists a, fi & [/in /ii+i] so that 9i{9c) = 6c + [a„] where 6c = 

6c{0,fi)- This is because when traverses [fin \ 6i{6c) traverses the interval 

{6i + [a„], 6r + [a„]). Let a be the value of a for fi. To solve for a fixed point we let 

^+ [aJ =6 + ai-uj(3^Hn¥ 

(4-16) 

z = hW"^ 

to obtain 
(4.17) 

Z = he'^ ^Q(a-[a„])^ 

From the first line we have 

(4.18) 1 + c sin ^ + E{6, jj) + a,.,a(/i) = e""'''(^-[""]). 

To solve (14.181) for 6, first we observe that 6c = 6c{0,fi) is a solution of (14.181) for 
z = 0. We then observe that 

Icos^el > K-\ 
This estimate follows from the fact that 6c is defined by 
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and F = e'^ '/3(a-[a„]) ^^^^ fUTTD . Applying the inverse value theorem to (I4.18p we 
obtain a solution 9 satisfying 

\9 - < Kfx. 

In summary we have obtained a fixed point (9, z) satisfying 

9^ 9c] z = be'^~'"(^-[^"]). 

To prove that {9, z) is an attracting fixed point, we compute the eigenvalues. The 
eigen-equation for DJ^ is 

- Tr{DT)\ + det(DJ^) = 0. 

From (I4.12P we have 
(4.19) 

det(DJ^) = al3-^W^'"-^^ « 1 



dz 



where for the first inequality we use 

d9i 





d^9^ 




8^9, 


( 


+ 




89^ 


898z 



<Ki\9-9,\ + \z\) 

at {9, z) with 

K = max 

Note that, on the domain the maximum is taken, F > i. The rest of 04.191] are 
obvious. It follows from (14.191) that both eigenvalues of DJ^ are close to 0. □ 

Our next Theorem is about the existence of non-degenerate transversal homoclinic 
tangency. 

Theorem 3 (Homoclinic tangency). Let the assumptions be identical to that of 

TheoremUi Then for every n > given, there exists fi G [/in ■* , /i^+i] , the corresponding 
value for a we denote as a, such that 

(i) Tg, has a saddle fixed point, which we denote as g(a), so that iy"(g(a)) fl 
W^*(g(a)) contains a point of non- degenerate tangency. 

(a) Let g(a) be the continuous extension of q{a) for a sufficiently close to a. Then 
as a passes through a, iy"(g(a)) crosses W^{q{a)) at the tangential intersection point 
of (i) with a relative speed > | with respect to a in 9-direction. 

Proof: Our plan of proof is as follows. We know that J-'a induces a horseshoe 
of infinitely many symbols in V \Vf, creating many saddle fixed points. Pick one 
and denote it as q. We prove that q is continuously extended over the /i interval 
[fin \ fJ'n+i]^ which we denote as g(a). Let iy"(g(a)) be the unstable and iy(g(a)) be 
the stable manifold of q'(a). We prove that iy"(g(a)) fl Vf has a horizontal segment 
traversing Vf, which we denote as £"(a). We also prove that H^*(g(a)) has a vertical 
segment fully extended in V, which we denote as £*(a). Observe that jFa(£"(a)) has 
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a sharp quadratic turn, and as fi varies from fin^ to fin\i, it moves from one side of 
V to the other, transversally crossing i^{sL). See Fig. 7. 




Fig. 7 Transversal homochnic tangency. 



Detailed proof for Theorem [3] is long and include some tedious computations. A 
complete proof is included in the Appendices. □ 

The following is a direct consequence of Theorem [31 

Corollary 4.1. Let the assumptions be identical to that of TheoremUl Then inside of 
every parameter interval [fin\ fJ-n+i], there is a set of parameters of positive Lebesgue 
measure, such that the homoclinic tangle associated with these parameters admits 
strange attractors with SRB measures. 

Proof: This follows from Theorem [3] applying |MV] and [BYj . both are based on 

[BC],toJ-a. □ 



4.3. Homoclinic tangles and observable chaos. Let J-'a be as in (14. 3 p and fi, A 
be as in (14. 9p . Q represents all solutions that stay close to i in forward times, and 
A represents all solutions that stay close to i in both the forward and the backward 
times. In this subsection we study numerically the structures of Q and A. 

We start with a concept of observability in numerical simulations. We say that 
a homoclinic tangle is observable in phase space if Q has positive Lebesgue mea- 
sure. Otherwise we say that this homoclinic tangle is not observable. We only expect 
observable homoclinic tangles to show up in numerical simulations. For maps with 
parameters, there is also an issue of observability in parameter space: a sub-collection 
of maps is observable only if it is from a parameter set of positive Lebesgue measure. 
See |WQk] for more detailed discussions on observable dynamical scenarios in numer- 
ical simulations. 

To numerically study homoclinic tangles through JF^, we drop the error terms in 
(14. 3 p and re-write kz as z. We obtain from (I4.3p - (I4.6I) a family of 2D maps in the 
form of 

= e + a - d ln(l + csm9 + z) 
zi=b[l+csm9 + z]\ 



'We thank Marcelo Viana for assuring us that, with Theorem [3l 'MVj directly apphes. 
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where a,b,c,d,'y are parameters, a G 5*^ b « 1, c > 1, d G M and 7 > 1. From 
Theorems and Corollary 14.11 we would expect at least three dynamical scenarios 
that are observable in parameter space. They are as follows. 

(i) For parameters of Theorem [T|, A for JF^ is a uniformly hyperbolic invariant set, 
and Q is the stable manifold for A inside of S~. Both Q and A are Lebesgue measure 
zero sets. The corresponding homoclinic tangle for these parameters is therefore not 
observable in phase space. 

(ii) For parameters of Theorem [21 Q contains an open neighborhood of a periodic 
sink so the associated homoclinic tangle is observable. Plots of individual orbits from 
Q would lead us to periodic sinks in A. 

(iii) For parameters of Corollary 14. ![ we expect strange attractors with SRB mea- 
sures to show up as an observable phenomenon in numerical simulations. 

In Figs 8-10, we plot Q and A for the maps defined in (14.201) with various choices of 
parameters that reflect the scenarios (i)-(iii) above respectively. Fig. 8 is for scenario 
(i), with a = 0.2, b = 0.005, c = 3, d = 2 and 7 = ^/2. Fig. 8(a) is a plot of all points 
in V, the orbits of which remain inside of V after 3 iterations. Fig. 8(b) is for after 
6 iterations. Nothing is left in V after 15 iterations. 



1] ! 
U 

'0 1.570! i.1416 4.7124 6.2!!2 'O 1.570! i.1416 4.7124 6.2!!2 

0<e<2]i 0<e<2]i 

(a) (b) 

Fig. 8 Homoclinic tangles with no sinks nor observable chaos, 
(a = 0.2, b = 0.005, c = 3, c/ = 2 and 7 = ^2) 

Fig. 9 is for scenario (ii) with a = 2. The values for b, c, d, 7 are kept the same 
as in Fig. 8. In this case A contains a periodic sink with a relatively large basin. 
Fig. 9(a) is for Q. All orbits initiated from Q quickly converge to an attracting 
periodic orbit. In Fig. 9(b) we depict 6^ v.s. k for one orbit from Q. Picture for Zk 
v.s. k is similar. Only one periodic sink shows up for A in numerical simulations. 
The horseshoe associated with the singularity of the logarithmal function (Smale's 
horseshoe), though exists inside of A, does not show up because the set it attracts is 
a set of zero Lesbegue measure. 
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Fig. 9 Homoclinic tangles with an attracting periodic solution, 
(a = 2, 6 = 0.005, c = 3, = 2 and 7 = ^2) 



Fig. 10 is for scenario (iii) with a = 1.5. The values for b, c, d and 7 are kept the 
same as before, fl is depicted in Fig. 10(a). In Fig. 10(b) we depicted again 9k v.s. 
k for one orbit. As k moves forward, 6k jumps randomly in a fixed range. These 
pictures represent a strange attractor with an SRB measure associated to transversal 
homoclinic tangency of a saddle periodic orbit of relatively high period. All orbits 
from Q in fact offer the same picture. 




Fig. 10 Tangles with an observable chaos, 
(a = 1.5, b = 0.005, c = 3, c/ = 2 and 7 = v^) 



We also performed systematic search over all combinations of parameters with b 
reasonably small. We persistently run into one of the three scenarios above. In the 
case of Fig. 8, however, sometimes it takes much longer for all points to be completely 
iterated out of V. This is particularly the case when d is small, and is more or less 
expected: as the overall strength of expansions around the horseshoe of Theorem [1] 
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gets weaker, the points in V tends to linger longer inside of V before been pushed 
out into U. 



4.4. Dynamical scenarios for periodically perturbed homoclinic solutions. 

In this paragraph we summarize all that have been obtain so far for equation (12. 2p in 
|WOj . |LW] and in this paper through the return map of Proposition 13. 3[ Again, we 
let Q{t) = sint and assume (HI) and (H2). The forcing parameters are inside of 



where » Rp » R^ » 1. Let W be the stable and IV" be the unstable 
manifold of {x, y) = (0, 0) in the extended phase space. Various dynamics scenarios 
for different parts of P are illustrated in Fig. 11. Since our purpose is to provide an 
overview, only descriptive statements are presented. Rigorous formulations and their 
proofs are either directly included in |W0] . |LW] and in this paper, or obtained by 
reasonable modifications of existing text. 

1. There is a surface S* in P (See Fig. 11), such that for all parameters under S*, 
n 7^ and we have homoclinic tangles for equation (12.21) . The dynamics of 
these homoclinic tangles are studied in Sect. 14.21 of this paper. In particular, there 
are open sets of parameters, such that the entire homnoclinic tangle is one uniformly 
hyperbolic horseshoe. There are also parameters for periodic sinks, and parameters 
for non-degenerate, transversal homoclinic tangency. 



2. For parameters over S* , fl IV" = 0. The return maps are again defined 
through (I4.3l) - (l4.7p . but for these parameters F > on E~ so JF^ are well-defined on 
S^. These maps have been studied systematically in |WY4] . We know that 

(a) There is a surface S above S* for which the following holds. For all parameters 
in between 5* and 5**, JF^ admit global attractors in S~ that are chaotic in the 
sense that they all contain a horseshoe (See |LW] ) . If the forcing frequency uj 



P={(a;,p,/i): a7G(0,i?J, pe{R;\Rp), /iG(0,i?;i)} 



p Q 





Fig. 11 Dynamical scenarios in parameter space. 
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is reasonably large, then there is a positive measure set of forcing parameters 
such that the strange attractors are rank one attractors of |WYlj and [WY2j 
with SRB measures (See |W0] ). 

(b) There is a surface Q (See Fig. 11) for which the following holds. For any give 
set of parameters on the left of Q, S~ is attracted globally to a simple closed 
curve, on which the induced map conjugates to a circle diffeomorphisms. If 
the rotation number of this circle diffeomorphism is rational, then there are 
saddles and attracting periodic solutions. If the rotation number is irrational, 
then the solutions are quasi-periodic. We also know as a fact that there 
are positive measure sets of parameters such that the corresponding rotation 
numbers are irrational and the corresponding solutions are quasi-periodic. 

(c) What happens between Q and 5* are as follows: as we move from the left to the 
right in the tu-direction, larger forcing frequency first deforms, then breaks the 
attracting invariant curve, inducing sinks and saddles. The unstable manifolds 
of these induced saddles would eventually fold in ^-direction, intersecting the 
stable manifolds to create strange attractors and rank one chaos. As we go 
down in the p-direction, and are pulled gradually closer. Reflected in 
the return maps of Proposition 13.31 is the growing relevance of the expansions 
associated with the singularity of the logarithmical function of (14.31) . 

5. HOMOCLINIC TANGLES FOR GENERAL FORCING FUNCTIONS 

In this section we let Q{t) be an arbitrary periodic function of period 27r. We 
explain how the different choices of the forcing function Q{t) affect the dynamics of 
the associated homoclinic tangles. 

Let A, S{ijj) and C{uj) be the same as before (See (14. ip ). A, S{uj) and C{ijj) are 
independent of Q{t). We assume (HI) and (H2)(i) for equation (12. 2p and replace 
(H2)(ii) by (H3) below. 

(H3) There exists a constant ^ > so that 



as \uo\ — i> +00. 

(H3) is stronger than (H2)(ii). It requires that the magnitude of the Fourier trans- 



decays exponentially as \uj\ oo. We note that there is no lack of known systems 
satisfying (H3). 

Expanding Q{t) in Fourier series we write 



^SHu) + CH7) 



formation R{ijj) of the function 



R{s) =e(s)e-^o^W'^" 



oo 



(5.1) 




n=l 
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If the mean value of a Q{t) is not zero, we give it to p. So there is no loss of generality 
in assuming (15.11) . Let us assume in addition that 

(5.2) cl + sl^ 0. 

Let Q(t) be as in (15. ip satisfying (15.21) and assume (HI), (H2)(i) and (H3) for 
equation (12.21) . Parameters uj,p,e,fi are specified as follows. First we fix (arbitrarily) 
an G [j^Ruj, Ru}]- We then fix a value of p such that 

3 < ■ ^cf + ■ VC2(^) + S-^ico) < 9. 

After that we fix e sufficiently small so that, 

\^(e)-Md)\ «pA 

where 4>l{6) is from (I3.19P and (j){6) is obtained by replacing —L~, with — oo, +oo 
respectively in (j)L{d)- We also make e sufficiently small so that 



2 < ^ ■ sjcl + si ■ ^CUu;) + SUu) < 10. 

p (<< e) is the only parameter we allow to vary. In what follows JF^ : is 
the return maps of Proposition 13.31 induced by equation (12.20 . Recall that p = \np. 

Theorem 4. Let Q(t) be as in Ii5. 1\) satisfying Ii5.2^) and assume (HI), (H2)(i) and 
(H3) for equation li2.S\) . Let the parameters uj,p,e be specified as in the above. Then 

(a) there are infinitely many disjoint open intervals of p, accumulating at p = 0, so 
that the corresponding homoclinic tangles of equation ^2.2\) are reduced to one single 
horseshoe of infinitely many symbols; 

(b) in between each of these parameter intervals, there are values of p so that the 
homoclinic tangles of equation ^2.2) contain stable periodic solutions; and 

(c) there are also parameters in between where the homoclinic tangles admit non- 
degenerate transversal homoclinic tangency. 

Proof: We argue that our previous proofs for Theorems [HIS] remain valid for the 
current setups. By assuming (113), (15. 2p and uj G [g^g-R^, which is >> 1, we 
make the first order term for Q{t) dominate in 0l(^)- Let us recall that 

(I)l{0) = j H(s) Q{e + UJS + uL-)e~ ^W'^-rfs 
is a critical element of the return map J^p in Proposition 13. 3[ By definition 

/oo 
M(s) Q{e + UJS + uoL-)e- ^W'^-rfs 
-oo 
oo 

= ^/c^ + s'^ ■ \/C'^{nuj) + S'^{nuj) ■ sm{n9 — nuoL^ — 6^) 

n=l 

where 9n are constants completely determined by Cn, Sn,C{nu), S{nu). We re-write 
J-'p of Proposition 13.31 following the steps of Sect. 14. H using z for X and denoting 
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h, Zi) = J-'p{9,z) for {0,z) G S . From Proposition 13.31 we have 



(5.3) 

where 

(5.4) 

(5.5) 
with 

(5.6) 
and 

(5.7) 



UJ 



p 

zi = h[¥{d,z,^P 



UJ 



a = - ln/i-1 + ^(L+ + + ^ ln(£(l + 0{e))P^ALp) 
¥{9, z,p) = l + c(sin e + ^{9)) + kz + E{9, p) + Oe,zM 



UJ) 



¥{9,p) 



c = {AlpY^ ■ + s\ ■ ^C\uj)^S^{, 
\^={A^p)-\P^Ptr\\^0{e)) 



{A^p)-\Pt)-\\ + Pl)0.,p(1) + {A^P)-\m - MO)) 



n=2 



cl + sj 



IC^{nuj) + S^{nuj) 



sin(n^ — nujL — 9„ 



By (H3) and the assumption that uj > -^Roj >> 1, $(6*) is an added error term, 
toward which our previous proofs of Theorems [TH are indifferent . □ 

From fl5.3p - fl5.7l) for JF^ we see that (I4.20p is a prototype of return maps for all Q{ujt) 
provided that uj » 1. If the forcing frequency is lower, then ^{9) in fl5.7p remains 
an important part of (f)L{9) in jFp. It is then possible to have a number of disjoint 
vertical strips for V, and an equal number of vertical strips for U. The images of each 
of the V-components again wrap around S~ infinitely many times in 9 direction. We 
might, however, have more turns for J-'{V), as shown in Fig. 12. 




Fig. 12 Infinitely wrapped horseshoe maps for Q{t) in general. 
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We finish by presenting two more numerical pictures for Q and A. Tliese are for 
tfie maps assuming tlie form of fl4.20p . but with sin 6' replaced by sin^ + sin 36'. Fig. 
13 is for the case of an attracting periodic sink, with a = 1, b = 0.005, c = l,d = 2. 
fl is depicted in Fig. 13(a). 9k v.s. k for one orbit from Q is depicted in Fig. 13(b). 
All orbits in Q is attracted to a periodic sink. 

Fig. 14 is for a strange attractor with an SRB measure. The values for b, c, d and 7 
are kept the same as in Fig. 13, but a is changed to 0.5. fl is depicted in Fig. 14(a), 
and 6k v.s. k for one orbit from Q is depicted in Fig 14(b). This orbit is attracted to 
an SRB measure. 




Fig. 13 Homoclinic tangle with a periodic sink. 
(a = l,b = 0.005, c = 1, d = 2 and 7 = v^) 




Fig. 14 Homoclinic tangles with observable chaos, 
(a = 0.5, b = 0.005, c = 1, ci = 2 and 7 = ^2) 
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Appendix A. Proof of Theorem [3] 



In this appendix we prove Theorem [3]0 Let a„ be the value of a at /i = fin "^ and 
[a„] = a„ — a„mo(i(27r). Let a(yu) be the value of a at /i G [Ain^/^i+i]- We divide the 
proof of this theorem into the following steps. 

Step 1. Solving for hyperbolic fixed points For /i G [/in \ /i^'^j^], let m be an integer 
> ?>uj(3~^ and gm(a) = (6',m, -Zm) be the solution of the equations 



UJ 



(Al) 



e + [sLn]+ 27rm = 6 + a(/i) - - lnF(6', z, n) 

P 



z = h[¥{e,z,iJ,)]7. 

9m is determined by 

(A.2) ¥i9m,zm,fi) = e-"/^('^('^)-['^"l-2--), 

and 

(A. 3) Zm = be'^"^"'''*^''-*"''*"^^^'^™-*. 

Claim A.l. gm(a) = {Om,Zm) is saddle fixed point. 
Proof of Claim \A.1[ Recall that 

o9 oz 

det(D^) = a(3'^W^'"-^^. 

oz 



Observe that, from (1A.2P and the assumption that m > So;/? ^, F( 
It follows that |ccos^m| > 1, and 

d9i 



Z-m: fJ') < inn- 



d9 



> 101. 



This implies 

\Tr{DJ^)\ > 100. 

Observe that we also have det(i5jF) << 1. Therefore we have two eigenvalues, one is 
close to and the other is with magnitude > 1. 

We also have 



Claim A.2. Form > 3iul3-\ 

d9r^ 



da 



< 



100' 



da 



< Kh. 



Proof of Claim [Aj.- Estimate for ^ follows directly from (]A.3I) . To estimate ^ 
we take derivative with respect to a on both-side of (]A.2p and use F < to obtain 
Iccosftml > 1. 



%inus Claim rOl b). which we prove in Appendix B. 
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Step 2. The stable and the unstable manifold for In the rest of this proof we let 
m be the smallest integer > 3Ljj3~^. Denote g(a) = Let 

V = {{9,z)e V, F>F(gloo„^(a),/i)} 

We obtain V from V by taking away two thin vertical strips at the vertical boundaries 
of V. Observe that, by definition, g(a) G V. We make e sufficiently small so that 
(1) the distance from g(a) to the vertical boundary of V is >> k; and (2) Km '■= 
IF(gioom) >> k. 

Denote the stable and the unstable manifold of g = g(a) as W^{q) and W^lq) 
respectively. The local stable and the local unstable manifold are denoted as 
and Wl^^{q). Let £"(g) be the connected branch of W^{q) in V \Vf that contains 
Wi^^{q), and ii{q) = J?-a(£"(g)). is a horizontal curve traversing Vf in 6'-direction: 
it is straight forward to verify that (1) i^{q) is a horizontal curve, (2) the image of 
is also a horizontal curve, and (3) the length of that image at least doubles the 
length of i'^{q) so it traverses Vf. Let z = w'^{9) be such that {9,w^{9)) G ^i(g). 



Claim A.3 

(a) 



We have on ii{q), 



(b) 



— I < bt ■ 

d^ 



Proof of Claim U.gl - 
have from fl4.12p 



(A.4) 



dw''{9i) 




d9i 





Denote = J^^. For {9,z) G r(g), let {9i,zi) = T{9,z). We 



.ld¥ dw^{e) 

dz de 



/I - .n-iid¥\ I , ,n-iid¥ dw^{e) 



(a) holds because the magnitude of the denominator in (]A.4p is > 1 for {9, z) G V\Vf. 



Remember that since t^{q) is horizontal we have 
derivative one more time to obtain 



de 



< 



100' 



To prove (b) we take 



(A.5) 



d^w''{9i] 
d^9i 



d_ 

de 



dw'^jei) 
dei 



dBi 

de 



Observe that is the denominator in (1A.4I) . the magnitude of which is > 1. Let 



M 



max 

(e,z)(^q 



d^w''{9) 



d9^ 



We have from flA.4l) and (lA.SP 



So 



and M < Kb < h^. 



d^zi 



<Kih + K2hM. 



d^9i 
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Let itiiq) be the segment of W^i^q) in V that contains Wf^^^q), and i^{q) = 

^{C-'Li{q))- C--i{q) is a fully extended vertical curve in V", which we represent by 
a function = w'^{z). 

Claim A. 4. We have on t_i{q), 



(a) 
(b) 



dw" (z) 



dz 
d'^w''{z) 



dz^ 



< ka ; and 



Proof of Claim\A^- Let {6, z) G ^(g) and denote = T'^iQ, z). We have 

from fl4^ 
(A.6) 



dd. 



dz 
dz—i 

dz 



dw'jz) 
dz 
1 



and it follows that 
(A.7) 
Let 

We have from flA.7p . 



because F > and 



dw^{z^i) 




a/J-ibF^^-'-if''";^") 




dz-i 









Ml = max 

{e,z)€etAi) 



dw^{z) 



dz 



dw^{z- 



dz-i 



< Ki\i + K2hMi 



It then follows that Mi < fTk < k 2 . 
To prove (b) we write 



(A.8) 



d w'^iZ-l) dz \ dz^i 



d ( dw^z-i) 



dz^i 



dz-i 
dz 



where '^^^j^^ is as in (1A.7P and 

dz_i 



(A9) 



Let 



dz a(3~ 



1 ( ^,^^u^o.,-^.id¥ dw^{z) ,,.-i^9¥\ 



Mo = max 



(e,2)e£li{g) 

We have from f[0) . flA8|) and flA9l) that 



< Kh{KihM2 + i^a) 
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from which we obtain M2 < Kh < ba . 







Step 3. Non-degenerate, transversal tangency Let £^ be a connected segment of 
£"(g) n Vf, and £^ be the vertical curve where ii{q), are as in Step 2. 

We use z = w^{9) to represent i'^ and 9 = w^{z) to represent if^. J^a(^a) traverses 
V in horizontal direction as /i runs through [/i^T^ /ii+i]- Consequently there exists 
ft G [fin \ fi^^\i], the corresponding value for a we denote as a, so that £| and J^a(^a) 
intersect tangentially at a point we denote as g = {9,z). Let (6*0,2:0) £ such 
that [9, z) = J-'a{9o, Zo). Our next claim implies that the tangential intersection of £| 
and J^a(^a) ^ is not degenerate. 



Claim A.5. For {9, z) G 1% let (^1, z^) = J^i 

(P9i 



z). Then at {9,z) = {9o,zo), we have 



dzf 



» 1. 



Proof of Claim \A.5[ From (14.121) we have 

d9, ii-cop-^m + up-^'-i'^ 



(A.IO) 



dzi a/3-ibF°/3-i-i|| + a/3-ibF"/3- 



dz de 



At the point of tangential intersection, we have 
unless 



d0i 



dzi 



< k2, which is not possible 



(A.ll) 



¥d9 



< h* 



from flA.lOl) . This is because ^"'^g^-' < b^ from Claim lAlST a). The effect of b in the 
denominator can not be possibly balanced if (lA.lip is false. 
For the estimate on second derivative we start from 



(A.12) 



where ^ is the denominator in 



de 



dzi 



d^9. 



f]~2 dzi 



To compute ^ ^^^^ derivative of 



the function on the right hand side of (lA.lOp with respect to 9. Applying the quotient 
rule we obtain a fraction, the bottom of which has a factor b^. On the top, we have 
a collection of finitely many terms, each of which is < Kh^^^, in magnitude except 
one in the form of 



(A.13) 



0-1^ d¥\ 



Remember that we have F > Km on V, and ^ > K ^ from flA.lip . We also have 



d_ 

d9 



'-"^ Wd9) 



de 
uj(3- 



F2 



fc^ + csin6') > 1. 
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Therefore, flA.lSp is the dominating term on top and we obtain 



dzf 

at q. <0 
To finish our proof of Theorem [21 we also need to prove that, as a varies, if^ and 
J^a{^^) move with different speed at the point of tangency. To make the dependency 
on parameter a exphcit, we write w'^ = w'^{6, a), = w^{z, a). Claim ETSl apphes to 
w^{6,sl) and Claim [A~4l applies to w'^{z,a). 



Claim A.6. Let 



2 . 



, a), zi{9o, a)) = J-'a{Oo, w^{9o, a)). Then at a = a we have 



(a) ||^^i(^o,a)| > f; and 

(b) \lw^{z,a)l<l^. 

Recall that q = {6, z) is the point of tangential intersection and (6*0, Zq) is such that 
J^k{0o,zo) = q. 

Proof of Claim IXT^ - In this prove we use dz, de and da to denote partial derivative 
with respect to 2;, 6* and a respectively. 

To prove (a) we let q = {9m, Zm) be the saddle fixed point and £^{q) and ^"(g) = 
^a(^"(g)) be as in Claim El For {9,z) e t{q) and (^0,^0) = ^a(^,^)- We have 
from 

00 = f{e,z,ai) = e + ai-'^\n¥{e,z,fi) 
(A.14) P 

zo=g{e,z,ai) = h[¥{e,z,fip, 

and in (lA.14p . zq = w^{9Q,a), z = w^{9,a) because both (^oj-^o) and {0,z) are on 
ii{q). We first invert the first equality in (]A.14p . obtaining 6 = 6{6Q,a.); then we put 
it into the second equality in (]A.14p to obtain zq = w^{9q, a). To estimate daW^{9o, a), 
we first let 

Ma= max |9aw"(^,a)| 
and obtain from the first equality in ( 1A.14I) . 

daf + dJ-daW^ 



\daeieo,a) 



<Ki+ K^Ma 



def + dj ■ dew- 

because \dof\ > 1 for (9, z) E V \ Vf and l^gw"! < from Claim lATsT a). From the 
second equality in flA.141) we have 

\daw\eo, a)| = \deg ■ daO + d^g ■ {dew^ ■ daO + daW^) + dag\ 

< KMKi + K^Ma) + K^h, 

from which it follows that 

(A.15) Ma < bi 

(a) now follows by taking da on 

01 (^0, a) = ^0 + a - uo(3^^ In F(^o, w;"(^o, a) , /i) 
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using flA.lSp . 



Proof of (b) is more sophisticated than that of (a). We need to study the stable 
manifold through the field of most contracted directions, a method originally intro- 
duced in |BC] and fully developed in [WYlj and |WY2j . A detailed proof is included 
in Appendix B. 

With Claim we know that, as a varies passing a, J^a(C) crosses £^ transversally. 
This finishes our proof of Theorem [3] owing that of Claim |X]6]^b). □ 

Appendix B. Proof of Claim [X?6](b) 

In order to produce the desired estimates in Claim lAlGT b). we need more precise 
controls on the stable manifold of the saddle fixed point q^- The main idea of our 
proof, that is, to approximate the stable manifold by using the integral curves of vector 
field defined by the most contracted directions of the Jacobi matrix, was originated 
from [BC], and was fully developed in [WYl] and |WY2j . Here we only need a specific 
version of the contents developed in the beginning part of Section 3 in [WY2] . 

B.l. Most contracted directions. In what follows Ui A U2 is the wedge product 
and {ui,U2) is the inner product for Ui,U2 G M^. 

Let M be a 2 X 2 matrix and assume M ^ cO where O is orthogonal and c G M. 
Then there is a unit vector e, uniquely defined up to a sign, that represents the 
most contracted direction of M, i.e. \Me\ < \Mu\ for all unit vectors u. From 
standard linear algebra, we know / = e"*" is the most expanded direction, meaning 
I Me"*" I > I Mm I for all unit vectors u, and Me ± Me"*". The numbers \Me\ and iMc"*"! 
are the singular values of M. 

Let u -L V he two unit vectors in M^. The following formulas are results of ele- 
mentary computations. First, we write down the squares of the singular values of 
M: 

(B.l) |Me|2 = - VB^ - AC) := A, \Mf\^ = + VB^ - AC) 

where 

(B.2) B=\Mu^^\Mv^, C=\MuKMvf. 

We write e = a^u + /5of, and solve for \Me\ = VX subject to + (3q = 1. There 
are two solutions (a vector and its negative): either e = ±f , or the solution with a 
positive u-component is given by 

(B.3) e = ^{au + f3v) 

with 

(B.4) a=\Mv^ -\ (3 = -{Mv, Mu) 

and 

(B.5) Z = y^a^ + /32. 
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From this we deduce that a solution for / is 

(B.6) f = ^(-pu + av). 

B.2. Stability of most contracted directions. In what follows we let qi = J^a{qo), 

M.-M^ i--r^f§ -r^lf 



We have for G V \ Vf, 

(B.7) 2 < \Ai\ < K, \Bi\ < K, \Ci\, |A| < Kb. 

Let M(") = DJ^^ilo)- = M„ ■ M„_i ■ ■ ■ Mi. Let the most contracted direction 

for M^") be e„ and the most expanded direction be /„. Denote the values of a,P 
and Z in (iRill and (!B3!1 for M^") as an, /3n and Observe that, assuming G 

(B.8) |M(-)/„| > L 

We have 

Lemma B.l. Let qq be such that qo, ■ ■ ■ ,qn ^ V \ Vf. Then for all 1 < i < n, 

(a) Id+i - dl < {Kh)\ |M«e„| < {KhY; 

(b) |aa(e,+i - e,)| < (i^b)\ |5aM«e„| < (irb)\ 

Proof: Let := jM^u A M^^^w]. We have 

(B.9) Ai = |det(M(*))| < (Kb)\ 

It then follows from |M«ei||MW/i| = A^ and flRSl) . 

(B.IO) |M«ei| < (i^b)'. 

We substitute m = Cj, f = fi and M = M(^+^) into dEl for e^+i and ([Rl]) for /i+i. 
By using ( IB. II) for M'^^^-^^/j+i, we have 

(B.ll) |M(*+i)/,| = iM^^+^Vml ± OiiKhY). 

from flB:2|) . (iRQl) and (iRTOj) . We also have 

(B.12) ^ ^ |M(^+i)/,|l 

We now prove Lemma [B. 1( a). Using u = Ci and v = fi, we have, from fIB.Sp . 

(B.13) Ci+i -ei = - — ^ — Ci + f3i+ifi 

To estimate |ej+i — ej|, we need to obtain a suitable upper bound for and lower 

bounds for \ai+i\ and Zi+i. We have from fiRil) . flBTOl) and flRT2l) . 

(B.14) lA+il < |M(^+i)e,||M(^+i)/,| < 
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and ~ ^j+i- These estimates together with Zj+i > 1 tell us 

The second assertion follows easily from 

|M«e„| < |M«(e„ - e„_i)| + ■ ■ ■ + |M«(e,+i - e,)| + |M«e,| < (irb)\ 
This finished our proof for Lemma [B. 1( a). 

To prove Lemma [B. 1( b) we start with 
Sublemma B.l. |9aei|, \dafi\ < Ki for some Ki. 

Proof: Let u = (0, 1)^, v = {l,Of and use flRS]) for d and flRGl) for /i. We have 
Zi > a > |Mif p — Kb > 1. Differentiating (1B.3P and (]B.6P gives the desired result. 



In the rest of this proof, d = d^- Our plan of proof for Lemma IB. 1( b) is as follows: 
For A; = 1, 2, ■ ■ ■ , we assume for alH < 

(*) |(9ej|, \dfi\ < 2Ki where Ki is as in Sublemma lB.il 

and prove for all i < k: 

(A) |9(M«/,)| < K\ |5(M«e,)| < {Khy- 

(B) |9(e.+i-e.)|, \dU\+i-fd\<iKhr- 

Observe that for i = 1, (*) is given by Sublemma lB.li It is easy to see that (B) above 
implies (*) with i = k + l, namely \dfk+i\ < \d{fk+i - /fc)| + ■ ■ ■ + |9(/2 -/i)| + \dfi\. 
From (B), we have — fi)\ < {Khy, and from Sublemma IB.lt we have \dfi\ < 

Ki. Hence \dfk+i\ < Kh + Ki, which, for b sufficiently small, is < 2Ki. The 
computation for e^+i is identical. 

Proof that (*) =^ (A): First we prove the estimate for d{M^^^ fi). Writing 

i 

9(M«/,) = ^ M, ■ ■ ■ {dM,) ■ ■ ■ Ml/, + M»9/,, 
i=i 

we obtain easily 

i 

|5(M«/,)| <Y,\M^■■■ {dMj) ■ ■ ■ Mi/,| + ||M«|| \df,\ < iK^ + K\2K{). 
i=i 

This estimate is used to estimate c}(M(*^ej). Write 9(MWej) = (/) + (//) where 
(/) is its component in the direction of M^^^ fi and (//) is its component orthogonal 
to M«/i. Recall that d{M^^ a, M^^ fi) = 0. We have 



|(//)| |M«/,| = |9(M«e,) A M«/,| < |9(M«e, A M«/^)l + iM^^^e, A 9(M«/^)I- 
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The first term in the last line is < {Khy, noting that we have established \dei\, \dfi\ < 
2Ki] the second term is < {Khy ■ K^. This completes the proof of (A). <) 

To prove (B), we first compute some quantities associated with the next iterate. 
Substitute u = Ci.v = fi,M = M(^+^) in flRB- flRGl) . The following is a straightfor- 
ward computation. 

Sublemma B.2. Assume (*) and (A). Then for all i < k: 

(a) |9A,+i| < (irb)2(^+i); 

(b) < (i^b)^v^; 



(a) \da,+,\,\dZ,+i\<K^y^. 

Proof that (*), (A) =^ (B): We work with e^; the computation for /j is similar. 
From (23) we have d{ei+i - e,) = (///) + {IV) + {V) where 



!(///)! = |^(em-e.)az,+i| < ^^/^ ■ (i^b)' < {Khy- 

\{m\ = < J-(|5A+ii + iA+ip/.i)<(irbr; 

m\ = \^d( e.)| « {Khy. 

To estimate (III), we have used Sublemma IB. 2( c) and part (a) of Lemma iB.li To 
estimate (IV), we have used Sublemma IB. 2( b). (*) and < (■^)*- The estimate 

for (V) is easy. 

This completes the proof of Lemma IB. 1( b) . □ 

We also need to control the speed of change for the most contracted directions in 
V \ Vf. Let qo{s, a) be a curve in V\Vf parameterized by a parameter s and assume 
that 

lko(s,a)||c2 < K. 

Let M^^\s) = DJ^^iloi^^ ^))) e„(s) be the most contracted direction for M(")(s). 

Lemma B.2. Let Qq be such that qo, ■ ■ ■ ,qn ^ V \ Vf. Then for all 1 < i < n, 

(a) \ds{ei+i{s) - ei{s)y < {Khy, |9,M«(s)e„(s)| < {Khy; and 

(b) |5A(e^+i(s) - e,(s))| < {Khy, |aAM»(s)e„(s)| < {Khy. 

Proof: The proof for Lemma [B. 2( a) is identical to that of Lemma [B. 1( b). It suffices 
to regard all d as ds instead of d^- The estimate for the second derivatives is proved 
by a similar argument. Here we skip the details. □ 

B.3. Temporary stable curves and the stable manifold. In the rest of this 
proof we let rj = bio and denote = {{0,z) E A : \z\ < rj}. We view e„ as a 
vector field, defined where it makes sense, and let 7„(s) be the integral curve to e„ 
with 7„(0) = go- 
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Lemma B.3. Let qq = be the saddle fixed point of Theorem\^ and 7n(s) he the 
integral curve to satisfying 7^(0) = go -A-q- Then, for all n> 
(^) \^ai.Q) " •^Kq'o)! < for all q = 7n(s) and all i < n; 

(b) 7n(s) is a fully extended vertical curve in {V \ V}) fl Ar^; 

(c) |7n+l(s) -7„(s)|, |(9a7„+l(s) -9a7„(s)| < hTo . 

Proof: Lemma [B.3( a) follows directly from |M*^'^e„| < (i^b)* for alH < n (Lemma 
IB.l( a)). Denote = JF^. Let Bq be the ball of radius 2t] centered at go- is 
well defined on Bq, and substituting u = (0, 1)^, v = (1, 0)"^ into ( IB. 31) we obtain 
s(ei) > Kh^^. Let 71 = 7i(s) be the integral curve to Ci defined for s G {—2ri,2ri) 
with 7i(0) = go. 

To construct 72, let Bi be the ?7^-neighborhood of 71. For ^ E Bi, let ^' be a point 
in 71 with le - e'l < V'- Then - ^(go)| < \H0 ' + ' < 

Kri^ + Khrj < Krj^. This ensures that 62 is defined on all of Bi. Let 72 be the 
integral curve to 62 with 72(0) = go- We verify that 72 is defined on {—2ri,2ri) and 
runs alongside 71. More precisely, let t G [0, 1] and 

g(t, s) =7i(s)+t(72(s)-7i(s)). 

We have 

1^(72(5) -7i(s))| < 162(72(5)) - 61(72(5))! + |ei (72(5)) - ei(7i(s))| 

< |e2 - ei| + |9fei||72(s) - 7i(s)| 

< Kh + K\-f2{s) - -fi{s)\. 

Here we use \dtei\ < K. By Gronwall's inequality, |72(s) — 7i(s)| < Kb|s|e^l'^l, which 
is << for |s| < 2ri. This ensures that 72 remains in Bi and hence is well defined 
for all s G (—2//, 2//). 

In general, we inductively construct 7, by letting -Bj-i be the r^^-neighborhood 
of 7i_i in S. Then for all ^ e - J^^(go)| = - go| < Kt]^ for 

k < i. Thus e.j is well defined. Integrating and arguing as above, we obtain 7^ with 
l7i(-s) - 7i-i(s)| < {Khy^'^lsl « rj' for all s with |s| < 2?]. 

To estimate the derivative with respect to a, we let 

qit, s) = 7„(s) + t(7„+i(s) - 7„(s)). 

We have 

|^<9a(7n+l(s) - 7n(s))| < 1 5a(e„+i (7„+i (s) ) - 6^(7^+1 (s) )) | 

+ |<9a(e„(7„+i(s)) - e„(7„(s)))| 
< \en+i - e„| + |9f9ae„||7n+i(s) - 7n(s)| < Kt]'^. 



From this the second item of Lemma IB. 3( b) follows. 
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□ 



110 



B.4. The proof of Claim [A. 61 (b). We are ready to prove Claim IX!6] (b). First we 
note that at the point of tangency, \z\ « bw so it suffices for us to consider A^^ in 
the place of A with 77 = bio. 

From Lemma IB. 3( c). we know that 7„ 700 uniformly as n 00, and from 
Lemma FB. 3 ( a) we know that 700 is the stable manifold of go = Qm-, which we write as 

(B.15) ^ = ^oo(s,a), 2 

We also have 

(B.16) l^a^ools, a) -9^^1(5, a)|, \d^eUs,s.) - dJi{s,s.)\ < 2bi 

from Lemma [B. 2( c). 

Write 7oo using 6 = w'^{z,a.). We ffist solve s = Soo{z,sl) from the second item in 
fIB.lSp . then substitute to the ffist item in flB.151) to obtain 

w''{z,a) = 6'oo(soo(-2, a),a). 

Differentiating on both side, we obtain 

d^w%z,ai) = 9,6'oo(s,a)9aSoo(2,a) + (9a6'oo(s, a) 

9s6'oo(s,a)^^24^^ + 9a6'oo(s,a) 



dgZoo ('5, a 
d^w{{z,a.) + 0{h 
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where 6 = a) is the equation for 71. To obtain the last estimate we use flB.16P 
and the fact that dsZi{s, a) > |. 

To prove Claim lAlGV b). it now suffices for us to conffim that 

d^w{{z,a.) < ^. 

To verify this estimate we observe 

(B.17) \d^l.wl{z,^)\ = \d^s-\e,)\ < Kb, 

where s(ei) is the slope for ei. This follows from a direct computation using flB.3p . 
From (lETTll . 

\d^wl{z,a)\ < |(9aWi(2„,a)| + Khr] 

where daWf{zm,Si) is the value of dawl{z,a.) at go = Qm- We now use Claim [Al2] for 
5aWi(^m,a). □ 
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